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We investigate the principal chiral model between two and four dimensions by means of a non 
perturbative Wilson-like renormalization group equation. We are thus able to follow the evolution 
of the effective coupling constants within this whole range of dimensions without having recourse 
to any kind of small parameter expansion. This allows us to identify its three dimensional critical 
physics and to solve the long-standing discrepancy between the different perturbative approaches 
that characterizes the class of models to which the principal chiral model belongs. 
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The field theoretical approach to critical phenomena 
has provided one of the greatest achievements of theoret- 
ical physics in the past thirty years EJ. This program 
has been successfully appHed to the SO{N)/SO{N — 
1) spin systems for which a coherent picture has 
emerged from the various perturbative (weak-coupling, 
low-temperature, l/A^V,|esj)ansions a, as well as non 
perturbative methods oETlP. However, for many sys- 
tems, such a satisfactory understanding has not yet been 
reached. This is actually the case of almost all mod- 
els whose order parameter is not simply a vector and 
whose symmetry breaking pattern departs from the usual 
SO{N) / SO{N — U scheme. Among them are frustrated 
antiferronramf ts el^tll , the dipole Jocked phase of super- 
fluid ^He E^By, superconductors Vjt3, 



Josephson junc- 
tion arrays O , electroweak phase transition llZl'Ea , smec- 
tic liquid crystal etc. For such systems, the differ- 
ent perturbative predictions conflict. On the one hand, 
the weak coupling analysis performed around four di- 
mensions on the Landau-Ginzburg- Wilson (LGW) model 
generally predicts a flrst orde.r uhaac transition due to the 
lack of a stable flxed point yillMEJ. On the other hand, 
the low temperature expansion performed around two 
dimensions by means of the corresponding Non-Linear 
Sigma (NLct) model predicts a flxed point in d = 2 + e 
dimensions-jSo-that a second order phase transition is 
expected OESIlHl . Finally, the results obtained from nu- 
merical and experimental studies in three dimensions 
vary from first to second order depending on the system. 
The disagreement between the perturbative approaches 
prevents one from safely extrapolating any of the per- 
turbative results in d = 3. Therefore, contrary to the 
SO{N)/SO{N ~1) case, the usual perturbative methods 
fail to provide the general picture of the critical physics of 
these systems. This problem clearly requires non pertur- 
bative approaches. In this letter, we study the principal 
chiral model, the archetype of such a situation, by means 
of a Wilson exact renormalization group approach, a very 
efficient method to go beyond perturbation theory. A co- 
herent scheme emerges, which might well be generic. The 
key point is that a non perturbative fixed point is found. 



the presence of which allows us to remove the discrepancy 
between the perturbative results in d = 2 and d = 4. For 
the principal chiral model, it leads to predict a first or- 
der phase transitieH in d = 3, in agreement with previous 
numerical works cJ . 

The principal chiral model pertains, aroQua others, to 
magnets with noBTcpUinear spin ordering 0'E2I as well as 
superfiuid '^He 113113. It describes the physics of systems 
whose order parameter consists in a triad of orthonor- 
mal three-component vectors R = (61,62,63) or, equiva- 
lently, in a rotation matrix. The lattice hamiltonian of 
the principal chiral model is given by that of 5*0(3) or 
0(3) rotation matrices, interacting ferromagnetically: 
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~JY.J2 <■< = -jJ2 Tr{'R^W} (1) 



<i,j> a—1 



<i.]> 



with J > 0, and deti?' = ±1 on each site. 

Note that for frustrated magnets and '^He the further 
constraint deti?* = -1-1 should or should not be added, 
depending on the underlying microscopic system. For the 
sake of simplicity, we discard this constraint in the fol- 
lowing. The Hamiltonian (|l|) is thus invariant under the 
50(3)00(3) group of left U £ 50(3) and right V e 0(3) 
global transformations: i?' UR^V. In the low temper- 
ature phase, the residual symmetry group consists in a 
diagonal 50(3) so that the symmetry breaking scheme is 
given by 50(3) ® 0(3) 50(3). The order parameter 
space of the model thus identifies with the 0(3) manifold. 

Beyond its physical interest, the principal chiral model 
provides - as a direct generalization of the 0{N) vector 
model to matrices - the simplest example lying in the 
class of models depicted above: its NLtO, [Sipdel analy- 
sis exhibits a fixed point in d = 2-f e mll3, so that 
a second order phase transition is expected, whereas 
the weak coupling expansion performed in d = 4 — e 
on the LGW model rather suggests a flrst order phase 
transition induced by fluctuations since no flxed point is 
identified Ej. In addition, the expansion predicts a 
continuous phasCr-transition everywhere between two and 
four dimensions in contradiction with the iV = 3 case 
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in d = 4 — e. Finally, the investigations by Monte Carlo 
simulations of the 0(3) model largely favor a propounced 
first order phase transition in three dimensions E3 . 

Let us now briefly come back to the different pertur- 
bative approaches, their domain of validity and draw- 
backs. As for the e = d — 2 expansion, the crucial fact 
is that the perturbative (3 function of a NLcr model asso- 
ciated with the symmetry breaking scheme G H only 
depends ouj-the local geometrical structure of the man- 
ifold G/H O. It is therefore insensitive to its topol- 
ogy. For the principal chiral model, the 0(3) mani- 
fold is locally isomorphic to the four dimensional sphere 
S'3 = SO{4:)/SO{3) whereas two manifolds globally dif- 
fer since no(0(3)) = ni(0(3)) = Z2, being free of 
topological defects. Therefore, the principal chiral model 
and the N p.4^ector model have the same perturbative 
/3 function E^^E^ although their low energy spectra differ 
by topological excitations. This raises serious doubts on 
the validity of the predictions based on the NLcr model 
approach, at least at sufficiently high temperature where 
these excitations are activated. On the other hand, the 
perturbative f3 function of the LGW model near d = 4 is 
not "geometrical" since it depends on the representation 
of G spanned by the order parameter and could be sensi- 
tive to the topological degrees of freedom. Let us however 
emphasize that there is no general consensus on the abil- 
ity of the perturbative treatment of the LGW model to 
deal with these topological degrees of freedom. In our 
case, both the LGW model and Monte Carlo simulations 
predict a first order phase transition, but this agreement 
should not be overestimated, since the absence of a fixed 
point near four dimensions is generally a poor informa- 
tion about the three dimensional physics. Indeed sev- 
eral systems, among which type-II superconductors, are 
suspected to undseap-a second order phase transition in 
three dimensions HO although no fixed point is found 
around four dimensions. For such mo de -I s ,, a non per- 
turbative approach in three dimensions BEHI has proven 
the unreliability of the e — A~ d expansion, extrapolated 
to e of order 1. The same situation is encountered for 
electroweak phase transition H . All these features show 
the necessity to go beyond perturbation theory. 

We perform here a non perturbative approach of this 
problem by means of an exact Wilson-like renormaliza- 
tion group equation— We use the concept of effective av- 
erage action Tk[<l>\ E3'E|1, which corresponds to a coarse- 
grained free energy functional and is obtained by inte- 
grating out all Fourier components of the classical field 
(j) with momenta q > k. The usual effective action - gen- 
erating one particle-irreducible correlation functions - is 
then recovered in the limit /c — *■ 0. The fc-dependence 
of Ffc is described by an exact evolution equatiop^ 
Legendre-transformed of the Polchinski equation 



SFfe 
dt 



4- 



dt 



(2) 



where t = Infc/A and A is an ultra-violet cut-off. The 
trace has to be understood as a momenta integral as well 



as a summation over internal indices. In Eq.(g), Rk is 
the effective infrared cut-off which suppresses the prop- 
agation of modes with— momenta q < k. A convenient 
cut-off is provided by : 
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(3) 



where is the field renormalization. In Eq.(|2|), F^.'^^' is 
the ea;acf/jeZ(i-(iepen(ient inverse propagator- i.e. the sec- 
ond derivative of the effective average action F^. This is 
what provides the non perturbative content of the equa- 
tion and its efficiency to recover critical-exponents of 
the three dimensional Heisenberg model Ej-eI, to inves- 
tigate the Kosterlilz^^hase transition O , the physics_of 
superconductors first order phase-transition 

low energy QCD H , Yang-Mills theory O , etc. 

Let us now deal with the principal chiral model. The 
exact evolution equation (|^) is a nonlinear functional dif- 
ferential equation, far too difficult to be solved exactly. 
We must truncate F^ in order to keep a manageable num- 
ber of coupling constants. We consider here a general 
S0{3) (8) 0(3) quartic effective action: 



Ffc 



1 A 2 

d'^X <! -Zk di_,<t)abd^^(j)ab + (/O - 3k) + /iT 



(4) 



where (j)ab is a real 3x3 matrix while p = Tr and 
r = i(Tr V0)^ - 5Tr(V0)^ are both 5-0(3) ® 0(3) in- 
variants. In Eq.(^ we have neglected all terms with more 
than two derivatives and derivative terms such as d^pd^p, 
assumed to be less relevant. For A > and p < 0, 
the minimum of the theory is realized by a configura- 
tion of the form f/)™" — \/l^Rab, Rab being a rotation 
matrix. The symmetry of the minimum is a diagonal 
50(3) group so that the symmetry breaking scheme is 
50(3)(g)0(3)/50(3). Note that for 0„f, = (/.™" one has: 
p — 3k and r = so that Eq.(^) corresponds to a quar- 
tic expansion around this minimum. The spectrum of 
the theory in the broken phase is obtained by studying 
small fiuctuations around 1^ consists in a triplet 

of Goldstone modes - corresponding to the three broken 
generators of the 0(3) group - and six massive modes re- 
lated to the fiuctuations of length and of relative angles 
of the axes of the matrix 0^;,: a quintuplet with mass 
mi = — 4/iK and a singlet with mass m2 = 2Xk . 

The coupling constants k, A and p, in the ordered 
phase, are defined by: 



Sp 



= 0; 



A = 6- 



6p^ 



St 



(5) 



Their fc-dependence is obtained by taking the i-derivative 
of these relations, by using Eq.(^ and the ansatz Eq.(p. 
The fiow for the field-renormalization Zk, rj = dlnZk/dt 
is computed along the same line (cf. ref. ]39[|). 
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The fixed points siiow up for dimensionless renormal- 
ized quantities: 



Kr = Zk "^K ; Xr = Z/^/ k 2 A ; /if 



k 2 ^ 



(6) 



In terms of tliese parameters, we liave obtained, after 
some work, tlie recursion equations in the ordered phase: 



dt 



dXr 

~dt 



-{d — 2 + ri)Kr — Vd 



= {d-4: + 2j])Xr - Vd 



2lf{0) + 2lf{mr,)- 



2Xll^{0) + 6Xll^{mr,)+ 



d^r 

~dF 



= (d - 4 + 2ri)fj,r + Vd A*r 



2 K 



4wd 
3d 



{Xr + 2/ir) Kr 2{Xr + 2/1^) Kr 

20/i^Krm2 2(0, fTlri) + 4A^Kr'Ti2.2 (0; "'■r2 ) 

(7) 



where v^^ = 2'^^^n'^/^r{d/2), m^^ — —ifirKr and nir^ — 
2Xr Kr- In Eq.(^ we have used the dimensionless "thresh- 
old functions" 



(8) 



^2,2(0,'^) 



dy y 



d/2 



d 



(1 + r + yr'Y 



dt [y{l + r)ny{l + r)+uj\ 



with y = q^/k^, r = 1/(6^' - 1) and s = -2ye?'/(e2' - if. 
In Eq.(^), acts only on the regulator r{y). These 

threshold functions encode the non perturbative aspect 
of the recursion equations (|^). As we now see, they also 
control the decoupling of the massive modes of the theory. 

The physics around two dimensions. Let us first show 
how our recursion equations allow to recover the predic- 
tions of the NLcT model around two dimensions. In the 
vicinity of d = 2, we expect the existence of a regime 
where the physics is controlled by Goldstone modes, i.e. 
where the massive fluctuations are frozen. This regime is 
defined through the limit of large mass: and — > 
oo. Using the dominant behaviour of m2^^(0,ct;) when 
uj ^ oo: m2^^(0,a;) ~ ~2uj^^, one gets: 



r] : 



SiTKr 



(9) 



On the other hand, the dominant contribution to the flow 
of Kr in this limit is given by: 



^^.(d-2 + 7j)Kr-2Vdim 

dt 



(10) 



since I'^i^j) oc 0;"^"+^^ when uj ^ 00. Using lf^'^{0) — —2 
and combining Eq.(^ and Eq.([lO|) one obtains the flow 
for Kr around two dimensions: 



dt 



1 

8^ 



(11) 



with e = d — 2. This is exactly the expected result pro- 
vided by the low temperature, expansion of the principal 
chiral model to first order IcjO. Note that the corre- 
spondance with the "temperature" T of the NLcr model 
is obtained by the substitution Kr — > 1/T. Equation ( pT| ) 
predicts a fixed point at a temperature T* = 1/k* = Stte 
of the SO{4:)/SO{3) universality class. We however 
emphasize that, if the existence of this fixed point is 
trustable in the immediate vicinity of two dimensions, 
where a Goldstone mode expansion is safe, it might be 
destabilized by non perturbative effects away from two 
dimensions where the critical temperature T* is finite. 

The physics around four dimensions. The weak- 
coupling perturbative expansion around four dimensions 
for the LGW action is recovered in the linear regime cor- 
responding to small masses: mr^,mr2 ^ 1. One then has 
to perform an expansion in powers of these masses. Using 
li{u;) ~ /^(O) - nwli+M for ^ « 1, and lt\Q) = -2 
we find 77 = 0{nl,Xl) and : 



dA,- 
"dT 

dflr 

~dr 



-tXr + 



87r2 



17 



xt 



40, 80 2 

— Xrllr + yMr 



(12) 



+ ^{-l2^ll+AXr^Mr) 



with e = 4 — d. This is the result found by means of a 
direct weak coupling expansion O . We recall that these 
equations admit no fixed point around d = 4. 

The physics between two and four dimensions. We 
have numerically analyzed the fiow equations (^ by vary- 
ing the dimension between d = 2 and d = 4. We have 
identified the different fixed points and studied their sta- 
bility. Note first that for /i = 0, Eq.(^ is the action of 
the 0(9)/0(8) vector model. As expected, we find the 
corresponding Wilson-Fisher fixed point, everywhere be- 
tween two and four dimensions. The situation is far more 
unusual for the principal chiral model. Around two di- 
mensions, we locate a stable fixed point C+ characterized 
by K* ~ l/87re (see Fig.l) and 77 ~ 3e with e = d — 2, 
which identifies with the standard S'0(4)/50(3) NLcr 
model fixed point. The coupling constants A* and /i* 
being finite at leading order in e, the masses to^ and 
diverge like 1/e, which justifies the Goldstone mode 
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treatment for this model. This situation remains quahta- 
tively unchanged until the dimension di = 8/3 is reached. 
At this particular value of d emerges a new and unstable 
fixed point C_ (Fig.l) whose coordinates behave like : 
K* ~ 1/e', A* ~ e'^ and /i* ~ est, where e' = d - 8/3, 



which corresponds to nir 



and 



0. This 



behaviour of the coupling constants as well as the ratio- 
nal value 8/3 can be obtained analytically from the flow 
equations (^. As the dimension d is increased, the two 
fixed points get closer, coalesce and finally disappear at 
d2 — 2.83. Above this dimension no fixed point is found. 
The fiow exhibits a runaway towards an instability re- 
gion, which can be interpreted as a first order behaviour. 
This result extends up to four dimen.sions, in agreement 
with the weak coupling expansion O . 

In conclusion, our approach, based on the concept of 
average action, allows us to overcome the diflticulties of 
perturbation theory and to obtain a coherent picture of 
the physics of the principal chiral model everywhere be- 
tween two and four dimensions. Qiiy results agree with all 
previous perturbative results oo^Eil and coincide with 
the numerical prediction .©f a first order phase transi- 
tion in three dimensions & . Previous studies of other 
models performed with the same method cil lead us to 
believe on the robustness of our results beyond our trun- 
cation Eq.(U). We conjecture that the appearence of a 
fixed point and its subsequent annihilation with the chi- 
ral fixed point is a generic phenomenon. This can be sub- 
stantiated by the observation that for generalizations of 
the principal chiral model (nx p matrices of real or com- 
plex fields, scalar electrodynamics, etc) perturbative stud- 
ies performed around four dimensionSj-e}dii bi | t the same 
kind of annihilation of two fixed points yMBEiy. The 
simplest assumption is that for all those models the ana- 
log of the C+ fixed point reaches two dimensions and 
identifies with that of the NLcr model whereas the other 
fixed point C_ vanishes with an infinite value of k* in 
a non trivial dimension d > 2. Finally, the precise di- 
mension c?2 where the two fixed points annihilate should 
be model-dependent. Whether it is greater, smaller or 
slightly smaller than three will imply the order of the 
transition in = 3 to be second, first or almost second 
order. 

We thank J. Vidal for a careful reading of the 
manuscript. 
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FIG. 1. K* as a function of d for the fixed points C+ and 
The fixed point C_ exists only between di = 8/3 and 
d2 — 2.83 where it annihilates with the principal chiral fixed 
point C+. Above d2 no fixed point is found. 



